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1. Introduction. 

Quantum Chromo Dynamics (QCD) with massless quarks of Nf flavours has SU (Nf) x SU (Nf) 
chiral symmetry which is dynamically broken at low temperatures by the vacuum. Since the u and 
d quarks have almost degenerate mass which in turn is much smaller than the scale of QCD, and 
the s quark is only moderately heavy, approximate chiral symmetry is a reasonable assumption for 
our world, i.e., QCD with 2+1 flavours of dynamical quarks. Rich phenomenological studies have 
demonstrated the utility, and consistency of such an assumption. Pions and Kaons are thus widely 
regarded as the Goldstone bosons resulting from the dynamical symmetry breaking by vacuum. 
It is thus a natural consequence that the finite temperature transition in our world is also widely 
accepted to be governed by chiral symmetry. It is therefore desirable to have the chiral invariance 
in any formulation which aims at studying the thermodynamics of QCD. Let us cite below two 
important physics aspects where it is even necessary to have exact chiral invariance. 

Lattice QCD has clearly been the method of choice for reliable non-perturbative studies in 
general, and finite temperature/density investigations in particular. As is well known, the Fermion 
doubling problem makes it difficult to demand the same chiral invariance as in continuum QCD. 
The original solution to this problem, namely, the Wilson Fermions, breaks all chiral symmetries. 
Since the staggered Fermions do have an exact chiral symmetry on the lattice, albeit at the cost of 
breaking of flavour and spin symmetries, they have dominated the area of nonzero temperatures 
and densities. 

As presented [jl|] in Lattice 2006 by one of us, the hadronic screening lengths illustrate their 
deficiency in the pionic screening length whereas the Overlap Fermions, with exact chiral, flavour 
and spin symmetry for any arbitrary lattice spacing, appeal - to do better. Advocated as means 
to explore the large scale composition of QGP, various numerical investigations with staggered 
quarks found that except the pion (and the sigma) screening length, all others could be understood 
as multiples of the appropriate number of free quarks (or antiquarks). The pionic correlator on the 
other hand showed nontrivial structure, and the pionic screening length approached the ideal gas 
value only in the continuum limit. The simulations with overlap quarks yielded a pleasant surprise 
on both counts and the corresponding pion screening length was close to ideal gas value for even 
small temporal lattices. 

Another fundamental aspect of QCD is the existence and location of the critical point in the 
T-ji. B plane, where /x B is the chemical potential for baryon number. Based on symmetries and a 
variety models, the QCD phase diagram is expected to have a critical point for two light and one 
moderately heavy quark. Again chiral symmetry plays a crucial role in this : the transition at fig = 
should be second order for two massless quarks, which turns into a cross over for light quarks. A 
line of first order phase transitions at finite density ought to terminate in a critical point. 

2. Ginsparg-Wilson Relation and ji ^ 0. 

Exact chiral invariance for a lattice Fermion operator D is assured if it satisfies the Ginsparg- 
Wilson relation [Q] : {/5,D} = aDy^D. In particular, the chiral transformations [^] 8\j/ = 0:75(1 — 
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fD)y and 8\j/ = 0!1/a(1 — jD)y%, leave the action S = £ X v \if(x)D X) y\if(y) invariant: 

55 = a £ % [ 7s D + Dy 5 - ^D 7s D - ^D 7s D] ^, = (2.1) 

The Overlap Fermions, and the Domain Wall Fermions in the limit of large fifth dimension satisfy 
this relation. Thus these Fermions with exact chiral invariance on the lattice are ideal for studies 
of the QCD phase diagram. One needs to introduce chemical potential in their known actions to 
do so. For the staggered quarks, this was done by first finding an expression for the conserved 
number, N, on the lattice, and then adding /J.N term to the action. It turned out [Q] that this lead to 
/1 2 -dependent divergences in the continuum limit even for the free case, which were removed by 
further modification of the action. Since the non-locality of the Overlap or Domain Wall Fermions 
makes the construction of a conserved charge difficult, it was proposed [§] that the Wilson Dirac 
operator Dw in the actions for these Fermions be modified using the same prescription as in the 
staggered case. This amounted to multiplying the link variable in the positive (negative) time 
direction by K(afi) = exp(a i u) {L{ap,) = exp(— a/i)). As the y$D-w{apL) is no longer Hermitian, 
the D OY = 1 + Y5Sgn(y 5 D w ) definition of the overlap operator necessitated an extension of the sign 
function: For complex X = (x + iy) eigenvalue, sgn(X) = sgn(x). Using this operator, it was shown 
[|]] showed numerically that no pL 2 -divergences exist in the free case (U = I). We demonstrated [f7|] 
the absence of the divergence in the free case analytically for all K,L such that K(a[i) -L(apl) = 1. 

We also claim [^] though that the chiral invariance is lost for nonzero jj,. This is easy to see by 
varying the quark fields again by the same infinitesimal chiral transformation as above : 

SS(an) = aY J %[y5D{apL)+D{apL)y 5 -^D{0)y 5 D{apL)-^D{apL)y 5 D{0)]^ y , (2.2) 

x,y 

whereas the extended sign function definition of the Dirac operator merely ensures 

y 5 D(an) + D(an)y 5 - a D(aii)y 5 D(an) = . (2.3) 

This is clearly not sufficient to make SS = 0. This is true for both Overlap and Domain Wall 
Fermions and for any K,L. 

A direct consequence is that the much desired exact chiral symmetry on lattice is lost for any 
/I 7^ 0, real or imaginary. Thus a pL -dependent mass will be acquired by even massless quarks 
in the interacting theory. The behaviour of chiral condensate as a function of pL will therefore be 
necessarily smoothened, wiping out any chiral transition that may be present. Of course, depending 
on how strong the transition is, it will begin to show up for small enough a or large enough Nj- 
How large a computational effort that may mean is a priori not clear. Recall that in the Taylor 
expansion method [j|] to incorporate the effects of nonzero jJ,, all the coefficients are evaluated at 
/i = 0. These will not suffer from any such lack of chiral invariance but the series in p, will be 
smooth, and will always exhibit convergence for any p.. 

One may be tempted to modify the chiral transformation itself for nonzero /I by demanding 
8y = a 7 5(l — |D (a ju ))!//■ and 8\j/ = ay/(l — |D(a^)) 75 . Clearly, 8S(an) = in that case. This 
is, however, not permissible since ysD{a}i), i.e, the generator of the transformation is not Hermi- 
tian. Moreover, a symmetry transformations should not depend on the "external" tunable parameter 
jj.. Recall that the chemical potential is introduced for charges Ni with [H,Nj] = 0. At least, the 
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symmetry should therefore not change as }X does. The most damaging practical consequence of 
such a modification of the chiral transformation is that the restoration of chiral symmetry due to an 
increase in }i cannot be addressed at all. For T / with }i = 0, the symmetry group remains the 
same at each T, allowing a change in the chiral order parameter (yy) (am = 0, T) to be interpreted 
as a change in the vacuum, i.e, restoration of the dynamically broken symmetry. With }i -dependent 
chiral transformations, the symmetry groups are different at each jx, with no obvious relation be- 
tween the respective chiral condensates. It thus appears much more reasonable not to alter the 
transformations but to look for a better way to include the chemical potential in the Overlap Dirac 
operator. 



3. Free Overlap and Domain Wall Fermions: Analytical Results. 




Figure 1: Contours in the complex o-plane for fi = (left) and T = O(right) cases. 

We investigated thermodynamics of free overlap and domain wall Fermions both analytically 
and numerically to demonstrate that i) the negative mass parameter < M < 2 is irrelevant in the 
continuum limit, and ii) the absence of jx 2 -divergences for general K and L. Using our numerical 
results, we also obtained an optimal range for M in order to obtain small deviations from the 
continuum limit on coarse lattices. 

The energy density, pressure, quark number susceptibility etc. can be obtained from In 3f = 
In det D ov by taking T, V, and ji or equivalently a/\ and a and a\l, partial derivatives. Here V = N^a 3 
and T = 1/(^7*04). Since the Dirac operator is diagonal in momentum space, its eigenvalues can 
be used to compute 2f: 



1 



h 2 = £h 



- [sgn (jh 2 + h 2 ^j h 5 ± i\fh 2 \/^h 1 + hl , with 
sinapj, j = 1,2 and 3, hi, = —a sm^a^p^jaA, 



!=1 



and h=. 



M- 



j 



■cos(apj)\ —a[l —cos(a4P4)]/a4 . 



(3.1) 



It is straightforward to show that e = 3P for all a and 04. Note that the free energy on finite 
volume is in general not equal to the pressure \P\, and thus will have 0(1/V) corrections spoiling 
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the equality above. Hiding the spatial momentum /^-dependence in terms of known functions [g] 
g, d and /, the energy density on an Af 3 x Nj lattice is 



ea 



N?Ni 



pj.n 



N?N T 



Pj,» 



£ (g + COS (O n ) + sjd + 2gCOS (0„ 



(3.2) 



'1 — cosa>„) 



sin (O n (g + cos (o n ) 



d + 2g cos ft),, (<i + 2g cos <o n ) (/ + sin 2 (O n ) 

where &)„(= £24/74) are the Matsubara frequencies. As in the case of the continuum, sum over 
(On can be carried out using the contour technique in the complex ft)-plane. The left panel in 
the Figure [T] displays the contour chosen for the /i = case. The crosses denote the Matsub- 



ara frequencies, the circles denote the physical poles in eq. (3.3) corresponding to the zero of 
(/ + sin 2 to) and the dashed lines correspond to the cuts. Evaluating the integrals, one obtains 



ea 



Pi 



[exp(Nr sinfT 1 a/7) + 1]~ ! + £3 + £4 , where / = sm 2 (apj), and 
the contributions from the contours at the top and bottom are denoted by £3 4 respectively. In the 
continuum limit, the above expression reduces to the ideal gas energy density ess for all M, with 
the cuts moving away to infinity faster than the contours at the to and bottom. 

For the case of T = but /J. / 0, the contour is displayed in the right panel of Figure [j]. 
Essentially the same treatment goes through as above, if one defines K(n) +L(ju) = 27? cosh 6 and 
K(n) — L(n) = 2i?sinh0 and substitutes in eq. ( ^3| ) Rsin((O n — id) for sinft)„ and similarly for 
cos (O n - After the contour integral one obtains 

£a 4 = _L£ 27rResF( J R,a))@^(a/i)-L(a^)-2 v / 7) +J F(R,(o)d(0-J F(\,(o)d(0 

(3.3) 

From the above expression, one notices that R = K(aji) -L{ajx) = 1 ensures cancellation of 
the last two terms. For R ^ 1, expanding R in powers of pi, on the other hand, the jii 2 -divergence 
is explicitly seen in the continuum limit. Note the higher order terms in that case violate the 
Fermi surface condition in the ©-function and contribute for all ji on the lattice; they vanish in the 
continuum limit. From the ©-function in the expression above, it is also clear that K and L should 
be such that K{apL) —L{a}i) = 2a pL + ff(a 3 ) in order that the usual Fermi surface condition is 
recovered. Furthermore, K(0) = 1 = L(0) is necessary for consistency with the zero density Dirac 
operator. We have also shown that the above derivations go through for the Domain Wall Fermions 
[Q] and for the general case of T 7^ and p, 7^ for both the Overlap ^ and Domain wall Fermions. 



4. Free Overlap and Domain Wall Fermions: Numerical Results. 

Numerical evaluation of the physical observables is a simple evaluation of sums over all al- 
lowed momenta apj = (2n/N s )rij with tij =0, N s — 1 and 0:4/74 = (2n4 + \)h/Nt with 724 = 
0, Nt — 1. Fixing a large £ = LT =N s /Nt for thermodynamic limit, we let Nj — > 00 to obtain the re- 
sults in the continuum limit. Figure g shows Ae(jl,T) /T 4 = [e(jl,T) -e{0, T)]/T A for n/T = 0.5 
for the Overlap (top panel) and the Domain Wall (bottom) Fermions for £ = 4. The £ = 5 results 
have been checked to be indistinguishable from these. Both sets of results show small, ~ 2-3 %, 
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Figure 2: pi -dependent contribution to the energy density for Overlap (top) and Domain Wall (bottom) 
Fermions for fi/T = 0.5. 



deviations from the continuum result (= 0.127 for pi/T = 0.5) for Nj > 12 for a range of M, being 
1.50 < M < 1.60 for Overlap Fermions and 1.40 < M < 1.50 for the Domain Wall Fermions. Note 
that the canonical M = 1 choice yields a smoother 1 /Nj dependence but larger deviations. 
Quark number susceptibility, defined as 



X y 



d 2 lndetD\ 
dH 2 ) 



T,V 



(4.1) 



can be worked out for the Overlap Fermions to be 



X 



N T N 3 a 2 ^ 



(/i 2 /i4 + h/[h^ cos(a£>4 — iaju)) u 
s 4 (s-h 5 ) 2 



+ 



s 2 (s-h 5 ) 



(4.2) 
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The u and v in the expression above are 



U = 2(S — rt5j(/l4- h«53 — )+s ( 



d/24 

da\l 



dh* 



da\L dafX 
dh 5 



and 



(4.3) 



(2h 4 + h +/i5Cos(a/?4 — iail)) — cos(a/?4 — iail) + 2/14/15 sinfa/74 — iapL) 

aajX 



with s 2 = h 2 + /j 2 and the hj as defined in eq.(3. 1) but with a^PA — ► {®P4 — iafX). Similarly one can 
write down an expression for it for the Domain Wall Fermions. Figure |5| displays the results for 
the Overlap Fermions (top panel) and the Domain Wall Fermions (bottom panel) as a function of 
y/Nj. The behaviour is each case is similar to the corresponding energy density case in Figure ||. 
In particular, the same range of M seems to be optimal for both Ae and ^(ju = 0). 



2.4 
2.2 

2 
1.8 
1.6 
1.4 
1.2 

1 

0.8 



M=1.00 
M=1.55 — * 
M=1.60 — e 
M=1.65 — & 




0.01 0.02 0.03 0.04 0.05 0.06 0.07 

1/N-r 2 



1.35 ■ 



1.20 ■ 




1.05 ■ 



0.90 ■ 



0.75 



Figure 3: Quark number susceptibility as a function of 1 /Nj for Overlap (top) and Domain Wall (bottom) 
Fermions for n = 0. 
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5. Summary. 

Exact chiral symmetry for the quark fields on the lattice, without any violation of other sym- 
metries, such as the flavour or spin, is important for several investigations of QCD thermodynam- 
ics. Indeed, the existence of the critical point in the T-}Xb phase diagram of QCD, and the non- 
perturbative determination of its location, depend on it. The currently popular choice of staggered 
Fermions may not be fully adequate for such studies. The Overlap and Domain wall Fermions are 
ideally suited but unfortunately lose their chiral invariance on introduction of chemical potential 
in the Bloch-Wettig method and its generalizations, as we showed above. We also proved analyt- 
ically that no 2 -divergence exists in the continuum limit for both the Overlap and Domain Wall 
Fermions for the Bloch- Wetting method and an associated general class of functions K{ji) and 
L(jj.) with K(/j.) -L(ju) = 1. Our numerical results corroborate these findings, and lead us to an op- 
timal range for the irrelevant parameter M. For the choice of 1.5 < M < 1.6 (1.4 < M < 1.5), both 
the quark number susceptibility at /J. = and the }i -dependent part of the energy density exhibited 
the smallest deviations from the ideal gas limit for Nj>12 for Overlap (Domain Wall) Fermions. 
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